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We present an efficient scheme which couples any desig- 
nated pair of spins in heteronuclear spin systems. The scheme 
is based on the existence of Hadamard matrices. For a system 
of n spins with pairwise coupling, the scheme concatenates 
cn intervals of system evolution and uses at most cn^ pulses 
where c « 1. Our results demonstrate that, in many systems, 
selective recoupling is possible with linear overhead, contrary 
to common speculation that exponential effort is always re- 
quired. 



I. INTRODUCTION 

In recent proposals to perform quantum computation 
in nuclear spin systems using nuclear magnetic resonance 
(NMR) techniques |]l]-Q, coupled logic operations are 
performed using spin-spin couplings that occur naturally 
in molecular systems. While this is straightforward for 
small systems with a few spins, generalization to complex 
molecular structures has been challenging. The compli- 
cation is caused by the many spin-spin couplings which 
occur along with the desired one. This fundamental task 
to turn off spurious evolution is so difficult that, coercing 
a complex system to do nothing - ceasing all evolu- 
tion - can be just as difficult as making it do something 
computationally useful. 

The task of turning off all couplings is known in the 
art of NMR as decoupling; doing this for all but a se- 
lect subset of couplings is known as selective recoupling. 
A common method to perform these tasks is to inter- 
rupt the free evolution by carefully chosen pulses. These 
pulses are single spin-1/2 (qubit) operations that trans- 
form the hamiltonian in the time between pulses in such 
a manner that unwanted evolutions in consecutive time 
intervals cancel out each other. 

Pulse sequences which perform selective recoupling are 
generally difficult to find for a large system. Each pulse 



simultaneously affects many coupling terms in the hamil- 
tonian. To turn off all but one of the coupling terms, 
these pulses have to satisfy many simultaneous require- 
ments. Ingenious sequences have been found for usual 
NMR applications but they do not address the 

problems relevant to quantum computation. In usual 
NMR applications, the structure of the spin systems is 
not known a-priori. Therefore, pulse sequences are de- 
signed to address all the spins together rather than indi- 
vidual spins. Quantum computation brings new require- 
ments, and initial efforts have been made to develop 
pulse sequences to satisfy these needs; however, to-date, 
schemes have necessitated resources (such as total num- 
ber of pulses applied) exponential in the number of spins 
being controlled. 

In this paper, we present an efficient scheme to per- 
form selective recoupling. In contrast to the situation 
with traditional NMR, this scheme addresses the prob- 
lem relevant to NMR quantum computation, in which 
the molecular structure is assumed to be well-known, and 
spins are individually addressible. The method is related 
to a class of well-known matrices called Hadamard ma- 
trices. We derive from any n x n Hadamard matrix a 
pulse scheme that decouples n spins using n time in- 
tervals and 0{n^) pulses. This decoupling scheme can 
easily be modified (i) to remove Zeeman evolution and 
(ii) to implement selective recoupling. This completes 
the construction for n spins whenever n x n Hadamard 
matrices exist. When n x n Hadamard matrices do not 
exist, a scheme for n spins can still be constructed using 
larger existing Hadamard matrices. In doing so, an ex- 
tra amount of effort is required, but this can be bounded 
using existence properties of Hadamard matrices. Al- 
together, the scheme requires cn time intervals and less 
than cn^ pulses, where c « 1 with upper bound c < 2. 
Our method applies whenever the spins couple pairwise 
and have very different Zeeman frequencies, such as in 
heteronuclear spin systems. 

The paper is structured as follows. In Section |l|, we 
review relevant concepts in NMR quantum computing 
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and re-state the problem precisely. In Section [II, we 
first motivate the construction of the decoupling scheme 
with examples, and then describe the general construc- 
tion related to Hadamard matrices. Important properties 
of Hadamard matrices are summarized. Modifications of 
the decoupling scheme to perform selective recoupling are 
described. We conclude with some general remarks and 
discussions of various properties and limitations of the 
scheme. 



II. NMR QUANTUM COMPUTING AND THE 
STATEMENT OF THE PROBLEM 



In this section, we describe the NMR system and de- 
scribe how a universal set of (non-fault tolerant) opera- 
tions [^ - [T^ , namely, the single qubit operations and the 
controUed-NOT gate, can be realized using basic NMR 
primitives. 

We shall consider a physical system which consists of 
a solution of identical molecules. Each molecule has n 
non-magnetically equivalent nuclear spins which serve as 
qubits. A static magnetic field is applied externally along 
the +z direction. This magnetic field splits the energy 
levels of the spin states aligned with and against it. This 
is described in the hamiltonian by the Zeeman terms, 
which, in the energy eigenbasis, are given by 
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where i is the spin index, LUi/2Tr is the Zeeman frequency 
for the i-th spin, and ai^^ is the Pauli matrix operating 
on the i-th spin. The convention ?i = 1 is used for the 
rest of the paper. The spins have very different Zeeman 
frequencies, a situation loosely termed as "heteronuclear" 
in this paper. 

Nuclear spins can interact via the dipolar cou- 
pling [0,^. This is given by the hamiltonian 
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where fy denotes the unit displacement vector from the i- 
th to the j-th spin, and gfj denotes the coupling constant 
between them. Spin-spin coupling can also be mediated 
by coupling to electrons. This indirect coupling has a 
tensor part and a scalar part. The tensor part is usually 
of the same form as Hd- The scalar part is given by the 
hamiltonian 
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However, in the presence of a strong external magnetic 
field, only the secular part (terms that commute with 
Hz) is important For a heteronuclear system, the 

resulting coupling becomes 
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independent of the original form of coupling. 

Single qubit operations are performed by applying 
pulsed radio frequency (RF) magnetic fields along some 
directions t) perpendicular to the static field. To address 
the i-th spin, the frequency of the RF field is chosen to 
approximate lOij^'K. When the w^'s are very different, 
very short pulses can be used, so that during the pulses, 
all other evolutions are negligible except for the rotation 
operator e~'i'^*'' '' where d is proportional to the pulse 
duration and the power. The Lie group of all single qubit 
operations can be generated by rotations about x and y. 
Our scheme uses only rotations of = tt along i, which 
implement (up to an irrelevant overall phase) on the 
spins being addressed. We denote this operation by X, 
superscripted by the spin index whenever appropriate. 

Coupled operations such as controUed-phase-shift or 
controUed-NOT acting on the i-th and the j-th spins can 
be performed given the primitive. 
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For instance, a controUed-NOT from the i-th to the j-th 
spin can be implemented by 
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The ultimate goal is to be able to efficiently realize 
arbitrary quantum operations on an n spin system with 
arbitrary couplings. In this paper, we consider a more 
limited objective, which can now be stated precisely, us- 
ing the definitions of Eq.(|l|), Eq.(P, and Eq.(||): 

Given a heteronuclear system of n spins with 
free evolution £"'(^2+^'^)*, controlled using 
typical RF pulses, how can 2Zij be imple- 
mented efficiently? 



HI. CONSTRUCTION OF THE SCHEME 

We will first construct a decoupling scheme to remove 
the entire coupling term Ji^ in the total evolution. The 
scheme is derived from Hadamard matrices, which will 
be reviewed. Methods to remove Jiz and to implement 
selective recoupling will be described afterwards. 



If the molecules tumble fast and isotropically, dipolar 
coupling and indirect tensor coupling will be averaged 
away; otherwise, the physics can be more complicated. 
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A. Construction of the decoupling scheme 

To construct the decoupling scheme, we only consider 
He in the evolution. Effects of Hz can be included later 
since all matrix exponents commute. 

To motivate the general construction, we analyze the 
simplest example of decoupling two spins. The evolu- 
tion operator for an arbitrary duration t is given by 
^ = g-jg^^t'^i"®'^?'. The sequence of events, tX'-^^tX'^^\ 
known as refocusing in NMR, will first couple and then 
decouple the spins. This can be described mathemati- 
cally as: 
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where 6 — 5^2^- is obtained using Taylor series 

expansion of the matrix exponents and using — /, 
and Eq. (|lO|) is obtained using anticommutivity of and 

The essential features of this decoupling procedure are: 

1. The pair of X^'^^ negates the sign of ci^-' in the evolu- 
tion between the pulses. 

2. The X pulses make the signs of the matrices of the 
two spins disagree for exactly half of the time. 

3. Since the coupling is bilinear in the CTz matrices, it is 
unchanged (negated) when the signs of the CTz matrices 
agree (disagree). The X pulses therefore negate the cou- 
pling for exactly half of the time. 

4. Since the matrix exponents commute, negating the 
coupling for exactly half of the time suffices for the evo- 
lution to be cancelled out. 

The most crucial point leading to decoupling is that, 
the signs of the matrices of the coupled spins, con- 
trolled by pairs of X pulses, disagree for half of the time. 

In general, we consider pulse sequences which concate- 
nate equal-time intervals and use X pulses to control the 
signs of the of each spin. The essential information 
on the signs can be represented by a "sign matrix" de- 
fined as follows. The "sign matrix" of a pulse scheme for 
n-spins with m time intervals is the n x m matrix with 
the (i,a) entry being the sign of cri*'' in the a-th time 
interval. These sign matrices have one-to-one correspon- 
dence with our restricted class of pulse sequences. We 
denote any sign matrix for n spins by Sn- For example, 
the sequence in Eq.(^ can be represented by the sign 
matrix 



The general construction of decoupling scheme is now 
reduced to finding sign matrices such that every two rows 
disagree in half of the entries. 

As a second example, we construct a decoupling 
scheme for four spins. We first find a correct sign matrix 
following the previous observations, and then derive the 
corresponding pulse sequence. For example, a possible 
sign matrix is given by 
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in which any two rows disagree in exactly two entries. S4 
can be converted to a pulse scheme by converting each 
column to a time interval before and after which X pulses 
are applied to spins (rows) given by — 's. No pulses are 
applied to spins (rows) with -|-'s. The resulting sequence. 



.(X(3)AWrA(3)x(4))(x(2)A(3)rA(2)x(3)) x 



(15) 



is the identity by construction and this can also be veri- 
fied directly. Note that Tic in r = e"*^"^* now denotes the 
sum of six possible coupling terms for four spins. Note 
also Eq. ( p^ is written in such a way that it corresponds 
visually to the sign matrix, though the evolutions are ac- 
tually in reverse time order relative to S^. However, such 
ordering is irrelevant for commuting evolutions. Since 



^2 = 



(13) 



= /, Eq.([T5|) can be simplified to 
T(x(3)A(*VA(4))(x(2)rX(3))(A('*)rA(2)x(4)) . (16) 

This simplified pulse sequence can also be obtained di- 
rectly from Eq. (^^ by converting columns to time inter- 
vals and inserting X^'-* between intervals whenever the 
«-th row changes sign or whenever a — sign reaches ei- 
ther end of the row. Pulse sequences for Eq.(p^ and 
Eq.(|l6|) are shown in Fig. 

The above scheme can be generalized to decouple n 
spins with m time intervals as follows: 

Construct the nx m sign matrix Sn, with en- 
tries -|- or — , such that any two rows disagree 
in exacly half of the entries. For each — sign 
in the «-th row and the a-th column, apply 
X'*^ before and after the a-th time interval. 

Because of the pulses, the sign of the cr^ matrix for each 
spin in each time interval is as given by the sign matrix. 
The az matrices of any two spins therefore have opposite 
signs for half of the time, so that their coupling is negated 
for exactly half of the time, and the evolution is always 
cancelled. 

For n spins, n x m sign matrices which correspond to 
decoupling schemes do not necessarily exist for arbitrary 
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FIG. 1. (a) Pulse sequence corresponding to Eq.(p^. From 
54, each "— " sign in the i-th row and a-th cohimn translates 
to two X pulses at cui before and after the a-th time inter- 
val, (b) Pulse sequence obtained from simplifying (a). This 
corresponds to Eq.(^6|), and can be constructed directly from 
S4 by translating each change of sign in the i-th row to an X 
pulse at oji. A "— " sign at the end of the row also gives rise 
to an X pulse at end of the last time interval. 



m, but they always exist for large and special values of 
m. A possible structure is: 



in which intervals are bifurcated when rows (spins) are 
added. Such bifurcation takes place whenever it is im- 
possible to add an extra row that is orthogonal to all the 
existing ones ( "depletion" ) . If such depletion occurs fre- 
quently, the sign matrix will have exponential number of 
columns, and decoupling will take exponential number of 
steps as n increases. The challenge is to find correct sign 
matrices with subexponential number of columns. This 
difhcult problem turns out to have solutions given by the 
Hadamard matrices, which will be decribed next. 



B. Hadamard matrices 

Hadamard matrices have applications in many areas 
such as the construction of designs, error correcting codes 
and Hadamard transformations JT^-pTt. 

A Hadamard matrix of order n, denoted by H[n), is 
an n X n matrix with entries ±1, such that 



H{n)H{n)^ = nl . 



(17) 



The rows are pairwise orthogonal, therefore any two rows 
agree in exactly half of the entries. Likewise columns are 
pariwise orthogonal. We abbreviate "±1" as "±". S2 
and S'4 in Eqs.(13) and ( p^ ) are simple examples of H{2) 
and i?(4). An example of H{12) is given by 
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The following is a list of useful facts about Hadamard 
matrices (details and proofs omitted): 

1. Equivalence Permutations or negations of rows or 
columns of Hadamard matrices leave the orthog- 
onality condition invariant. Two Hadamard ma- 
trices are equivalent if one can be transformed to 
the other by a series of such operations. Each 
Hadamard matrix is equivalent to a normalized one, 
which has only +'s in the first row and column. For 
instance, H{12) shown previously can be normal- 
ized by negating the 7-th row and column. 

2. Necessary conditions H{n) exists only for n — 1, 
n = 2 or n = mod 4. This is obvious if the matrix 
is normalized, and the columns are permuted so 
that the first three rows become: 

!-••• + + ■•• + + ■•• + + ••• ^ 

!-••• + + ■•• H ... - - ... - 



3. Hadamard's conjecture pq | H{n) exists for every 
n = mod 4. This famous conjecture is verified for 
aU n < 428. 



4. Sylvester's construction 19| If H{n) and H{m) 
exist, then H{nm) can be constructed as H{n) ® 
Him). In particular, HiT') can be constructed as 
iJ(2)®'', which is proportional to the matrix rep- 
resentation of the Hadamard transformation for r 
qubits. 

5. Paley's construction [po| Let q be an odd prime 
power. If g = 3 mod 4, then H{q + 1) exists; if 
g EE 1 mod 4, then H{2{q + 1)) exists. 

6. Numerical facts [|l4| For an arbitrary integer n, let 
n and n be the largest and smallest integers that 
satisfy n < n <n with known H{n) and H{n). We 
define the "gap", Sn, to be n — n (see Fig. For 
n < 1000, H{n) is known for every possible order 
except for 6 cases, and the maximum gap is 8. For 
n < 10000, H{n) is unknown for 192 possible orders 
and the maximum gap is 32. 
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FIG. 2. The gap 5„ between two existing orders of 
Hadamard matrices. 

The importance of the full connection to Hadamard 
matrices will become clear after we construct the scheme 
for an arbitrary number of spins in the next section. 

C. General Construction of Scheme 

In this section, a decoupling scheme for an arbitrary 
number of spins is constructed. Variations of the scheme 
to remove Zeeman evolution and to implement selective 
recoupling are also constructed. The requirements of the 
scheme are discussed. 

Recall that n is the minimum integer no smaller than 
n with known H(ji), and m is the number of columns in a 
sign matrix. For any variations of the scheme for n spins, 
77i„ is defined to be the minimum number of columns in 
any valid sign matrix and it represents the number of 
time intervals needed for the intended operation. 

To construct a decoupling scheme for n spins when 
H{n) does not necessarily exist, any H(n) with n — n 
rows omitted can be used as the sign matrix since subsets 
of rows of H{n) are still pairwise orthogonal. In other 
words, any nxn submatrix of H(n) is a valid sign matrix 
for decoupling n spins, and m„ = n. 

To remove both Hz and He, we use the fact that the 
Zeeman term for the i-th spin is linear in tri*'' , and negat- 
ing di*' for half of the time results in no net Zeeman evo- 
lution for the i-th spin. Therefore, Zeeman evolution for 
all spins can be removed if the sign matrix has identically 
zero row sum. Such a sign matrix can be constructed by 
starting with a normalized H(n) and excluding the first 
row of H(ri) in the sign matrix. Since a normalized H{n) 
has only -|-'s in the first row, all other rows have zero row 
sums by orthogonality. Such construction is possible un- 
less 71 = n, in which case construction should start with 
H{n + 1). Therefore, m„ = n ii n < n and m„ = n -I- 1 
otherwise. 

To implement selective recoupling between the i-th and 
the j-th spins, the sign matrix should have equal i-th and 
j-th rows but any other two rows should be orthogonal. 
The coupling term g^jal (E) cri, never changes sign and 
that coupling is implemented, while all other couplings 
are removed. The sign matrix can be obtained from a 
normalized H{n) by replacing the 1-st row with the j-th 
row, and replacing the j-th row with the i-th row. This 
scheme also removes Hz and r7i„ = n. To implement 



^ij, the duration of each interval t is chosen to satisfy 
g'^jUt = 7r/4. Note that the total time used to implement 
2Zij is the shortest possible, since the coupling is always 
"on" . 

The scheme requires m„ time intervals. Consequently, 
it requires at most nnin pulses, since XX — I and the 
X pulses are only used in pairs. The remaining ques- 
tion is: how does m„ depend on n? For simplicity, we 
consider n in place of m„. If Hadamard matrices exist 
and can be constructed for all orders, n = n. However, 
some Hadamard matrices are missing, either because no 
construction methods are known or they simply cannot 
exist. Due to missing Hadamard matrices, n = cn where 
c > 1. We argue in the following that the scheme is still 
very efiicient. First of all, we prove c < 2. For each n, 
there exists r such that 2^~^ < n < 2^ . Since H{2^) 
exists by Sylvester's construction, cn = n < 2^ < 2n. 
We now show that c is close to the ideal value 1 in most 
cases, due to the existence of Hadamard matrices of or- 
ders other than powers of 2. This is why the full connec- 
tion to Hadamard matrices is important. In Fig. H, c as a 
function of n is plotted for n < 10000. Within this tech- 
nologically relevant range of n, c deviates significantly 
from 1 only for a few exceptional values of n when n is 
small. One arrives at the same conclusion by considering 
the gap Sn = n — n > n — n, which bounds the extra num- 
ber of intervals needed in the scheme caused by missing 
Hadamard matrices. For n < 10000, Sn/n ■€i 1 except 
for the few exceptional values of n as a numerical fact. 
For completeness, we present arguments for c « 1 for ar- 
bitrarily large n in Appendix This is based on Paley's 
construction and the prime number theorem. Finally, if 
Hadamard's conjecture is proven, (5„ < 3 Vn. 




FIG. 3. Plots of c vs n, where cn = n = m„ is the minimun 
number of time intervals required to perform decoupling or 
selective recoupling for an n-spin system, c for n < 100 and 
101 >n< 10000 are plotted separately. 



IV. CONCLUSION 

We reduce the problem of decoupling and selective re- 
coupling in heteronuclear spin systems to finding sign 
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matrices which is further reduced to finding Hadamard 
matrices. While the most diSicuh task of constructing 
Hadamard matrices is not discussed in this paper, solu- 
tions already exist in the literature. Even more impor- 
tant is that the connection to Hadamard matrices results 
in very efficient schemes. 

Some properties of the scheme are as follows. First of 
all, the scheme is optimal in the following sense. The 
rows of Hadamard matrices and their negations form the 
codewords of first order Reed-MuUer codes, which are 
perfect codes ^6 17\. It follows that, for each Hadamard 
matrix, it is impossible to add an extra row which is or- 
thogonal to all the existing ones. Therefore, for a given n, 
nin = n is in fact the minimum number of time intervals 
necessary for decoupling or recoupling, if one restricts 
to the class of pulse sequences considered. Second, the 
scheme applies for arbitrary duration of the time inter- 
vals. This is a consequence of the commutivity of all the 
terms in the hamiltonian, which in turn comes from the 
large separations of the Zeeman frequencies compared to 
the coupling constants. Spin systems can be chosen to 
satisfy this condition. Finally, disjoint pairs of spins can 
couple in parallel. 

We outline possible simplifications of the scheme for 
systems with restricted range of coupling. For example, a 
linear spin system with n spins but only /c-nearest neigh- 
bor coupling can be decoupled by a scheme for k spins 
only. The i-th row of the nx k sign matrix can be chosen 
to be the r-th row of H{k), where i = r mod k. Selec- 
tive recoupling can be implemented using a decoupling 
scheme for k + 1 spins. The sign matrix is constructed as 
in decoupling using H{k + 1) but the rows for the spins 
to be coupled are chosen to be the k + 1-th row differ- 
ent from all existing rows pT[ |. This method involving 
periodic boundary conditions generalizes to other spatial 
structures. The size of the scheme depends on k and the 
exact spatial structure but not on n. 

The scheme has several limitations. First of all, it only 
applies to systems in which spins can be individually ad- 
dressed by short pulses and coupling has the simplified 
form given by Eq.(^. These conditions are essential to 
the simplicity of the scheme. They can all be satisfied 
if the Zeeman frequencies have large separations. Sec- 
ond, generalizations to include couplings of higher order 
than bilinear remain to be developed. Furthermore, in 
practice, RF pulses are inexact and have finite durations, 
leading to imperfect transformations and residual errors. 

The present discussion is only one example of a more 
general issue, that the naturally occuring hamiltonian in 
a system does not directly give rise to convenient quan- 
tum logic gates or other computations such as simulation 
of quantum systems . Efficient conversion of the given 
system hamiltonian to a useful form is necessary and is 
an important challenge for future research. 
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64-94. 

[24] H. Davenport, Multiplicative Number Theory, (Springer- 
Verlag, New York, 1967). 



APPENDIX A: UPPER BOUNDS FOR iV 

An argument for c w 1 for large n is pr esent ed using 
Paley's construction (mentioned in Section III B| ) , known 



results on primes in intervals and the prime number the- 
orem for arithmetic progressions. 

Let 7r(x) be the number of primes which satisfy 2 < 
p < X. For X < 67, a;/(logx — 1/2) < Tr{x) < x/{\ogx — 
3/2) 1^^. It follows that there exists a prime between n 
and n{l + e) for e > 2/logri. Applying Paley's construc- 
tion, H{p+1) or H{2{p+1)) exists depending on whether 
p = 3 mod 4 or p = 1 mod 4. Therefore, n < n{l + e) + 1 
or n < 2(n(l -I- e) -I- 1) respectively. 

The worse of the upper bounds n < 2{n{l + e) + 1) 
resulting from p = 1 mod 4 can be improved. Note that 
there are at least r primes between n and n(l -I- e)*". If 
the primes that equal 3 mod 4 and 1 mod 4 are randomly 
and uniformly distributed, the probability to find a prime 
which equals 3 mod 4 between n and n(l -|- e)*" is larger 
than 1 — 2^^. This assumption is true due to the prime 
number theorem for arithmetic progressions p^ . Let 
7r(a;, a, q) denotes the number of primes in the arithmetic 
progression {a, a + q, a + 2q, . . .} which satisfy 2 < p < x. 
It is known that 7r(a;, 3,4) w 71(2;, 1,4). Therefore, with 
probability larger than 1—2"'", n < n(l-f e)'"-t-l, implying 
c < (1 -|- e)'' -|- 1/n « 1 for large n. 
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